clear;
clc;
A=[16 2 3 13;
    5 11 10 8;
    9 7 6 12;
    4 14 15 1];
B=inv(A),
A*B
disp('--------------------------------------------------')
A=sym(A);
inv(A)
